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Abstract. Recently, surface magnetic field maps had been acquired for a small sam- 
ple of active M dwarfs, showing that fully convective stars (spectral types ~ M4 and 
later) host intense (~ kG), mainly axi-symmetrical poloidal fields. In particular, the 
rapidly rotating M dwarf V374 Peg (M4), believed to lie near the theoretical full convec- 
tion threshold, presents a stable magnetic topology on a time-scale of ~ 1 yr. The rapid 
rotation of V374 Peg (P = 0.44 days) along with its intense magnetic field point toward 
a magneto-centrifugally acceleration of a coronal wind. In this work, we aim at investi- 
gating the structure of the coronal magnetic field in the M dwarf V374 Peg by means of 
three-dimensional magnetohydrodynamical (MHD) numerical simulations of the coro- 
nal wind. For the first time, an observationally derived surface magnetic field map is 
implemented in MHD models of stellar winds for a low-mass star. We self-consistently 
take into consideration the interaction of the outflowing wind with the magnetic field 
and vice versa. Hence, from the interplay between magnetic forces and wind forces, 
we are able to determine the configuration of the magnetic field and the structure of 
the coronal winds. Our results enable us to evaluate the angular momentum loss of the 
rapidly rotating M dwarf V374 Peg. 



1. Introduction 



The rotational evolution of M dwarf (dM) stars can be inferred from observations of 
open clusters at different ages ( Irwin et al.l2006l: Scholz & Eisloffell2007l : lHartman et al. 



20091 : ICollier Cameron et aDl2009l : lMeibom et al.ll2009h . In voung (< 700 Myr) open 



clusters, dM stars still present high rotation rates, which suggests that angular momen- 
tum l osses at the early main-sequence phase are negligible for them (Ir win & Bouvieri 



2009). However, as the cluster ages (> 700 Myr), the number of rapidly rotating dM 



stars decreases, implying that there should exist a mechanism of angular momentum 



remo val that acts on time-scales of a few hundred million years (Sc holz & Eis loffel 



20071) . For solar-like main sequence stars, the magnetised stellar wind is believed to 
spin down the star by carrying away stellar angular momentum. It has been obser- 
vationally established that the a ngular velocity r ate Oq for solar-like stars varies as a 



function of age t as Q.q k t ^ 2 (Skumanichl l 19721) . However, it seems that the empirical 



ISkumanichl 's law is not valid for low-mass stars, suggesting that a solar-type wind (i.e., 
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with low velocities and mass-loss rates) cannot reproduce the rotational evolution of 
fully-convective stars. 

The existence of hot coronae, rapid rotation, and high levels of magnetic activity 
in dM stars suggests the presence of winds with an enhanced mass loss as compared 
to the solar wind. However, the low-density, optically thin winds of these stars pre- 
vents the observation of traditional mass-loss signatures, such as P Cygni profiles. The 
still unobserved high mass-loss rates from dM stars could be able to disperse debris 
discs, explaining why discs around dM stars older than > 10 Myr are scarcely found 
(IPlavchan et al.ll2005l) . Estimates of mass-loss rates from dM stars vary considerably. 
It has been suggested that the coronal winds of dM stars, despite of being very tenu- 
ous, possess mass-loss rates (M) that can co nsiderably exceed the solar value (M - 
2 x 1( T 14 Mr, yr 1 ) by factors of 10 to 10 4 (iMullan et all 1 19921: iBadalvan & Livshits 
19921 : iLim & Whitelll996l : Ivan den Oord & Dovld 1 19971 : 1 War serin & Drake! 1200 lh . al- 



though IWood et all (120011) claim an upper limit of M < 4 x 10~ 15 M Q yr 1 for Proxima 
Centauri (dM5.5e), 5 times below the value of the solar wind mass-loss rate. 

In this work, we investigate the coronal wind of a specific fully-conve ctive dM star, 
V374 Peg, for which obs erved surface magnetic maps have been acquired dDonati et all 
20061 : iMorin et all 12008). For this, we use three-dimensional (3D) magnetohydrody- 



namics (MHD) simul ations based on our previous m odels developed for solar-like stars 
d Vidotto et al.ll2009bh and weak-lined T Tauri stars dVidotto et al.ll2009al l2010bh . For 
the first time, an observationally derived surface magnetic field map is implemented 
in MHD models of stellar winds for a low-mass star. V374 Peg is a suitable case for 
modelling as a first step, because its surface magnetic distribution is close to potential, 
which implies that the adopted boundary conditions match the observed map closely. 
We self-consistently take into consideration the interaction of the outflowing wind with 
the magnetic field and vice-versa. Hence, from the interplay between magnetic forces 
and wind forces, we are able to determine the configuration of the mag netic field and 



the stru cture of its coronal wind. More details of this work can be found in lVidotto et al. 
(I2010ah . 



2. The Numerical Model and Results 

V374 Peg has mass M t = 0.28 M Q , radius R t = 0.34 R Q and is rotating with negli- 
gible differential r otation (i.e., as a solid body) with a period of rotation Pq = 0.44 d 



(IMorin et al.1 12008). We consider that its axis of rotation lies in the z-direction. To per- 



form the simulations, w e use the 3D MHD numerical code BATS-R-US developed at 
University of Michigan ( Powell et al.lll999l) . which solves the ideal MHD equations. 



The simulations are initialised with a ID hydrodynamical wind for a fully ionised 
plasma of hydrogen. Immersed in this wind we consider a magnetic field anchored on 
the stellar surface that has a geometry derived from extrapolations from observed sur- 
face magnetic maps using the potential-field source surface (PFSS) method (Figured^). 
The MHD solution is evolved in time from the initial magnetic field configuration to a 
fully self-consistent solution (Figure [Tp). We do not adopt fixed topologies for either 
the magnetic field or for the wind, as both the wind and magnetic field lines are allowed 
to interact with each other. 

In the PFSS model, the stellar wind plasma is not included directly, but its ef- 
fects on the magnetic field (and vice-versa) are incorporated through the inclusion of 
the source surface. Such a surface, for instance, alters the number of open magnetic 



On the Angular Momentum Loss of Low-Mass Stars: V374 Peg 



3 




Figure 1 . (a) Initial configuration of magnetic field lines for the extrapolation of 
the surface map using the potential-field source surface (PFSS) technique, where the 
grey circumference represents the position of the source surface, (b) Final configu- 
ration of the magnetic field lines after the self-consistent interaction with the stellar 
wind. 



field lines, through where a stellar wind could escape. The source surface (grey sphere 
in Figure [Ik) is chosen to lie at rss = 5 where beyond that, the magnetic field is 
considered to be purely radial. PFSS methods are usually criticised because their basic 
assumptions (the magnetic field is a potential field and the source surface is spheri- 
cal) may not always be met. However, the advantage of the PFSS method over the 
MHD models relies on its simplicity: it is simpler to implement and requires much 
less computer resources than MHD models. In our model, we use the magnetic field 
configuration derived by the PFSS method as initial condition and boundary condition 
at the surface of the star. We note that the surface of the star that occupies co-latitudes 
> 120° is never in view as the star rotates and so the magnetic field there can not be 
reconstructed reliably. 

Our simulations require a set of input parameters for the wind. Unfortunately, 
some of them are poorly constrained by observations. For V374 Peg, the magnetic 
field is the bet ter-constrained param eter. We have, therefore, implemented in our pre- 
vious models dVidotto et alj|2010bh surface magnetic maps derived by data acquired 
in 2005 Aug (iDonati et al.ll200fj|) . These observations show that V374 Peg hosts an in- 
tense, mainly axi-symmetrical dipolar magnetic field, with maximum intensity of about 
1660 G, i.e., 3 orders of magnitude larger than the Sun. 

The wind temperature and density are less constrained for V374 Peg. We, there- 
fore, adopt values representative of dM stars. dM stars are believed to host coronae 
with a high-te mperature plasma ~ 10 7 K in conjunction with a low-temperature one 2 
3 x 10 6 K dSchmitt et alJll990l : IcTiampapa et all 1 19961) . In our models, we adopt a 



temperature at the base of the coronal wind of Tq = 2 x 10 6 K or 10 7 K. These coronal 
temperatures are about the same order of magnitude as the solar coronal temperatures 
ofl.56xl0 6 K. 

Coronal densities inferred from X-r ay observations of dM stars suggest densities 
ranging from 10 10 cnT 3 to 5 x 10 12 cnT 3 dNess et al.ll2002l . l2004) . Therefore, we adopt, 



at the base of the coronal wind, densities in the range 10 10 - 10 12 cm . Compared to 
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the solar coronal density of about ~ 2 x 10 8 cm -3 , coronal densities inferred for dM 
stars are about 2-4 orders of magnitude larger than for the solar corona. 

The density, along with the magnetic field, are key parameters i n defining the 
magne tic field configuration of the stellar wind and its velocity profile dVidotto et al. 



2009b. a). Together, they define the plasma-/?, defined by the ratio of thermal to mag- 
netic energy densities. Therefore, at the base of the coronal wind of V374 Peg, 

fh = S^TT « 2 - 5 x l^ioTfi. CD 
^/(8tt) 

where the index "0" means the variable is evaluated at the base of the coronal wind, 
mo = 7i /(10 10 cm" 3 ) and T 6 = T /(2 x 10 6 K). For mp = 1 and Tf, = 1 Bp is about 



5 orders of magnitude smaller than for the solar wind dPneuman & Kopplll97lh . This 
implies that the winds of dM stars are highly magnetised and, therefore, are expected 
to differ from solar- type winds. 

Table [Qpresents the parameters adopted for the set of simulations we performed. 



Table 1. Adopted parameters for the simulations. The columns are, respectively: 
the case name, the density no and temperature To at the base of the coronal wind 
(r = R t ), the plasma-/? parameter evaluated at R* [Eq.[T|, the mass loss rate M, the 
angular momentum loss rate J, and the time-scale for rotational braking t. 



Case 


no 


To 


A> 


M 


j 


T 




[cm" 3 ] 


[MK] 




[10- 11 Moyr 1 ] 


[10 33 erg s- 1 ] 


[Myr] 


IMap 


10 10 


2 


2.52 x 10" 5 


4.2 


3.4 


180 


2Map 


10 11 


2 


2.52 x 10- 4 


14 


7.6 


84 


3Map 


10 12 


2 


2.52 x 10- 3 


50 


32 


17 


4Map 


10 11 


10 


1.26 x 10- 3 


26 


9.1 


48 



We were able to find a MHD solution for the wind for all the simulations we ran, 
showing that it is possible to develop coronal wind models with a realistic distribution 
of magnetic field. In general, MHD wind models are studied under the assumption of 
simplistic magnetic field configurations, especially when in pursuit of an analytical so- 
lution. Therefore, the study of a magnetised coronal wind where an observed magnetic 
field distribution is considered has long been awaited. Our work also sheds some light 
on the yet unobserved winds from dM stars. 

By comparing cases where only the base coronal density no was varied, we found 
that the poloidal velocity of the wind scales approximately as 

2 1 

u oc — for a given Bo- (2) 
no 



This qualitatively agrees with previous results ( Vidotto et all l2009bh . where it was 



found that an increase in the density leads to winds with lower velocities. Simula- 
tions presented here are in a very low-/? regime, where magnetic effects completely 
override thermal and kinematic effects of the wind. Therefore, Eq. © should be treated 
with caution, as under different B regimes (for example, when it approaches Bq ~ 1), it 
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becomes invalid. Figure^ presents the scaled wind velocity profile u p for cases IMap, 
2Map, and 3Map, 

1/2 fl\ 

u p - u p n{ 2 , (3) 
where n\2 = «o/(10 12 cm" 3 ). We note that the wind terminal velocity is Uco « [1300 

1 /2 1 

- 2100]ni 2 1/i km s , where the range of velocities refers to different wind latitudes 
(low-wind velocity near the equator, high-wind velocity around the poles). Because 
the magnetic field in the lower hemisphere of the star is not reliably reconstructed (co- 
latitudes > 120° of the surface of the star are not observed), a high-velocity wind 
develops there. Although this feature is local and does not affect the remaining parts of 
the grid other than radially away from the stellar surface, it is an artifact of our method 
and should not be taken into consideration (e.g., see the accumulation of magnetic field 
lines near the low-hemisphere of the star in Figures |2^ and[2j)). 

Equation Q also implies that the mass-loss rate of the wind (M oc pu r ) should 
scale as 

1 /2 

M oc pu r oc nJ , (4) 

which means that, despite the fact that the wind velocity of case 3Map is 10 times 
smaller than case IMap [Eq. ©], its mass-loss rate is one order of magnitude larger 
than for case IMap [Eq. ©]. This has implications on the efficiency of angular mo- 
mentum loss, as will be shown later. The mass-loss rates for cases IMap, 2Map, and 
3Map are M » 4 x \Q- W n^ M Q yr" 1 . 

Figure [2Jd shows the scaled poloidal velocity profile u p for case 4Map. This 
case considers a different temperature at the base of the corona (10 7 K as opposed 
to 2 x 10 6 K), and, because of that, has a larger /?o (Table [T]). For this case, we did 
not find an analytical expression relating velocity and temperature. The Alfven sur- 
face location and configuration of magnetic field lines are similar to the other dipolar 
cases, but the scaled wind velocity u p is smaller than the previous cases: ~ [850 - 

— 1/2 — I 

1600]?i 12 km s . The lower velocity observed in case 4Map happens because of its 
higher /3o- The mass-loss rate for case 4Map is M - 2.6 x 10~ 10 M Q yr -1 . 

Overall, our solutions differ considerably from the solar wind solution, where a 
low-velocity wind (terminal velocities of Uoo,q - 400 - 800 km s -1 ) with low mass- 
loss rate (M© - 2 x 1 0~ 14 M yr -1 ) is found . We note that, based on more simplistic 
wind models, such as Web er & Davisl ([1967), in the fast magnetic rotator limit, wind 



terminal velocities of 3320 km s are expected for a wind mass-loss rate of about 
10~ n Moyr 1 . 



3. Discussion: Angular Momentum Evolution 

Observations of the rotation evolution of dM stars in open clusters at different ages 
provide a way to constrain the time-scale t for the angular-mom entum loss. It has 
been suggested that r ~ 200 Myr or, mostly likely 400 - 800 Myr, dScholz & Eisloffel 
120071) for dM stars. Angular momentum of the star is carried away by the stellar wind. 
Because in our simulations there is no axi-symmetry, the torque j on the star has x, y 
and z components. Here, we are interested only on the z-component, as it is the only one 
responsible for the rotational braking (because the angular velocity of the star points in 
the z-direction). The z-component of the angular momentum carried by the wind is 
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Figure 2. Meridional cut of scaled poloidal wind velocity u p — u p / yjnn profile 

(a) for cases where a magnetic surface map was used IMap, 2Map, and 3Map, and 

(b) 4Map. Black lines represent the magnetic field configuration, and white line is 
the location of the Alfven surface. 



dMestel & Sellevlll970h 



J = 



+ J [r x (az x r)] z pV • MS A , 

Js A 



Jv A 



az x r x p(V + az x r)dV A 



(5) 



where V = u-o-zxr is the velocity vector in the frame rotating with angular velocity az, 
z is the unit vector that points in the z-direction, S a is the Alfven surface that delimits 
the volume Va, and n is the normal unit vector to the Alfven surface. The first term on 
the right of Eq. © does not contribute to the z-component torque and is therefore null. 
The second term disappears in the case of a spherical Alfven surface, but it is non-null 
in the cases where a surface magnetic map is considered and it becomes relatively more 
important for the cases with larger adopted fl<j. The third term is the dominant term in 
our simulations. 

We can estimate the time-scale for rotational braking as r = J/j, where J is the 
angular momentum of the star. If we assume a spherical star with a uniform density, 
then / - 2/5M t RlQ.o and the time-scale is 

where Pq = 2n/Q.o is the rotational period of the star. For V374 Peg, this results in 

6.45xl0 35 A)r 
t * : Myr. (7) 



Because / depends on the mass flux crossing a given surface, i.e., on the mass-loss 
rate of the wind M, from Eq. (0]), we have a rough scaling relation between j and M for 
cases IMap, 2Map, and 3Map 

joc Mcc n J /2 , (8) 
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which implies in a time-scale [Eq. (0] for rotational braking that scales as 

t cx n Q m . (9) 

— 1/2 

For cases IMap, 2Map, and 3Map, r 18n 12 M yr, well below the estimated solar 



spin-down time t© 7 Gyr (We ber & Davrl l967). 



Table [T] presents the mass and angular momentum loss rates, and the time-scale 
for rotational braking calculated for all simulations, where we verify the approximate 
scaling given by Eqs. (§]), ([8]>, and (©. Comparing to the observationally derived rota- 
tional braking time-scales of a couple of hundreds of Myr for dM stars is open clusters 
(IScholz & Eisloffelll2007h . we tend to rule out cases with larger coronal base densities 



(i.e., riQ > 10 11 cm" 3 ). According to this comparison, the most plausible wind density 
is the one used for models IMap. Such a density is also able to reproduce typical emis- 
sion measures of dM stars (EM ~ 10 51 cm -3 ) and comparatively (with the remaining 
cases) smaller mass-loss rates and higher wind velocities. Ultimately, when the star 
ages, the stellar rotation brakes, reducing the stellar surface magnetic field intensity, 
and therefore the wind velocity. 

With the inclusion of an observed distribution of surface magnetic field, we head 
towards a more realistic modelling of magnetised coronal winds. Never the less, our 
model presents limitations, such as the neglect of a detailed energy balance. Instead, 
we consider a polytropic relation between pressure and density parametrised through 
y in the derivation of the energy equation of the wind. Once the magnetic field dis- 
tribution is set, the thermal pressure adjusts itself in order to provi de a distribution of 
heatin g/cooling that is able to support the MHD solution obtained (lLow & Tsinganos 



119861) . If the wind of V374 Peg is able to cool down, e.g., by radiative cooling not con 
sidered in our models, the terminal velocities of the wind could be considerably smaller. 
Depending on where in the wind energy deposition (or removal) occurs, the wind ve- 
locity may change, without affecting the mass-loss rates. For instance, if a substantial 
cooling occurs above the Alfven surface, the velocity profile of the wind from that point 
outwards will be affected. As the information of what is happening above the Alfven 
point cannot be transmitted to the sub-Alfvenic region, the wind density and velocity 
profiles in the proximity of the star will not be changed, and consequently neither the 
stellar mass-loss/angular momentum-loss rates. 

Acknowledgments. The simulations presented here were performed at Columbia 
(NASA Ames). 
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